Near-optimal two-mode spin squeezing via feedback 
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We propose a feedback scheme for the production of two-mode spin squeezing. We determine a 
general expression for the optimal feedback, which is also applicable to the case of single-mode spin 
squeezing. The two- mode spin squeezed states obtained via this feedback are optimal for j — 1/2 
and are very close to optimal for j > 1/2. In addition, the master equation suggests a Hamiltonian 
that would produce two-mode spin squeezing without feedback, and is analogous to the two-axis 
countertwisting Hamiltonian in the single-mode case. 

PACS numbers: 03.65.Ud, 42.50.Dv, 32.80.-t, 03.67.-a 



I. INTRODUCTION 

Two-mode spin squeezed (TMSS) states are important 
to consider for a number of reasons. First, they may be 
used to demonstrate entanglement experimentally using 
only spin measurements [Q. Second, as they are entan- 
gled, they may be used as a resource for quantum infor- 
mation, for example, quantum teleportation quan- 
tum computing, superdense coding, etc. In addition, 
TMSS states are analogous to two-mode squeezed states 
for light , which have proven to be useful for quantum 
information purposes Q . 

There are several approaches to produce TMSS states. 
Recently they have been produced experimentally for 
samples of atoms by quantum nondemolition (QND) 
measurements of the spin components Q. In Ref. g| 
it was proposed to produce TMSS states of light by com- 
bining the two Einstein-Podolsky-Rosen (EPR) output 
modes of an optical parametric oscillator with coherent 
light at polarizing beam splitters. 

It is also possible to adapt proposals for generating 
single-mode spin squeezing to the two-mode case. For 
example, it has been proposed that single-mode spin 
squeezed states can be produced via an adiabatic ap- 
proach 1^. In this approach, one starts with a Hamilto- 
nian for which the ground state is easily achieved, and 
slowly varies the Hamiltonian towards one that has as its 
ground state the optimal spin squeezed state. In prin- 
ciple, this procedure can be adapted to obtain optimal 
TMSS states. A drawback to this method is that, in 
practice, the appropriate Hamiltonian may be difficult 
to realize. 

Another technique for generating TMSS states that 
may be generalized from the single-mode case is feed- 
back. Feedback has proven to be a powerful technique for 
state preparation ||, |^ and measurement jl0[ |l^ . It has 
been shown that it is possible to produce near-optimal 
spin squeezing in the single-mode case by continual mea- 
surement of the spin operator , and using feedback to 
bring (J^) close to zero 0. Here we adapt this method 
to the QND measurement approach discussed in Ref. 

The method for production of TMSS states in Ref. 
suffers the drawback that, although the variances in the 



sums of the spin components are small, their means are 
not close to zero. As is shown in Ref. the optimal 
TMSS states have zero means. We would therefore wish 
to obtain states with zero means in order for them to 
be as close as possible to the optimal states. This is 
analogous to the problem with the production of single- 
mode spin squeezing as considered by Thomsen, Mancini, 
and Wiseman (TMW) @. 



II. TWO-MODE SPIN SQUEEZING 

In Ref. 101, TMSS states (conditioned on the measure- 
ment record) are obtained via QND measurements of the 
spin of the atomic samples. In order to perform this 
QND measurement, an off-resonant pulse is transmitted 
through two cesium gas samples. The interaction Hamil- 
tonian may be expressed as M na 



Hint^ha I Sz{z,t)J^{z,t)dz, 



(1) 



where a is a coupling constant, Sz{z,t) is the z compo- 
nent of the instantaneous Stokes vector for the light, and 
Jz{z,t) is the z component of the continuous spin oper- 
ator. These variables are explained further in Appendix 

In the Schrodinger picture Jz{z,t) is independent oft, 
and the total Stokes vector S(z) is independent of z. The 
total time evolution operator is, therefore. 



[/ = exp ( -- / Hintdt 



exp y^ia J Sz{z,t)Jz{z)dz dt 

exp ^— za J SzJz{z)dz 
ejip{—iaSzJz)- 



(2) 



In Ref. the light is initially linearly polarized, so the 
light is in the maximally weighted Sx eigenstate. In 
addition, the total spin s is large, so we may apply a 
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SU(2)^HW(2) contraction. That is, we make the sub- 
stitutions 1^, |l^ 



Si, 



S^-]f,ib-,^), 



(3) 
(4) 

(5) 



where b^b is the number operator for photons subtracted 
from the x hnear polarization mode (and added to the 
y hnear polarization mode). We will assume that the 
probe beam is in a coherent spin state oriented in the x 
direction (i.e., x linearly polarized), and take the limit 
s ^ (X) so that this contraction is exact. 

Under this contraction, [Sy^Sz] — *■ is. This means 
that, under the Hcisenberg picture, the operator Sy 
transforms as 



QOUt Qin 



[iaSzJzj Sy] ~ S*'" H ttJz 



(6) 



where n = 2s is the total number of photons in the 
pulse. Now we consider two samples, where we will use 



the notation jj^^^ and J^'"' , where k € {x,y,z}, for the 
spin operators for samples 1 and 2, respectively, and 



t(2) 



J, 



(±) 



jI^^ ± J^.'"' . We may replace Jz with J^~^' , the 



r(2) 



(+) 



total z component of spin for the two samples, giving 



S?. 



(7) 



Thus a measurement of 5°"* 



gives a QND measurement 

of ji^''. In considering detection, it is more convenient 
to consider the Stokes vector for the light integrated over 
time interval 5t. The transformation of this component 
is 



5St 



a6n 



(8) 



where 5n is the total number of photons in time St. In 
Ref. TMSS states are obtained by performing QND 

measurements on Jy^^ as well. They achieve this by ap- 
plying a magnetic field in the x direction, in order to 
induce Larmor precession of the y and z components of 
the spin. This means that Eq. (^) becomes 

6S;^Ht) « SS'^it) + ^[Ji+) cos(m) + J(+) sin(f7t)], 

(9) 

where fl is the frequency of the precession. In this case 
measurement of JS*™* alternately gives QND measure- 
ments of ji^'' and J^^^ , thus giving reduced uncertainty 
in both Ji+^ and 

For the TMSS states considered in Ref. m there is the 



requirement that {J, 



(2)\ 



^ 1. Here we will 



instead take the convention that both {Jx^'') and (J; 
are large and close to the total spin j, and replace J, 



.(2)^ 
(+) 



with Jy . This is equivalent to a trivial rotation of co- 
ordinates for one of the modes and does not require a 
physical alteration to the experiment. If we now define 



ji^^t) = 4+^ cos(m) + j^-'> sm{nt), 

then Eq. (^ becomes 

aSn 
Y 



5s;-\t)^ssi-{t) + ^jfHt). 



(10) 



(11) 



In order for this to be accurate, we require Sn to be 
large, so that we may perform the SU(2)^IIW(2) con- 
traction. Under this contraction, 6Sy = \J 6s/2X, where 
X = b + h^ . Therefore the transformation in the quadra- 
ture X is 



X°"*(i) ^ X'"(<) + a^/5^J'p\t). 



(12) 



As the initial state is equivalent to the vacuum state [un- 
der the SU(2)— >HW(2) contraction], it has a variance of 
1 in its measured value. The same is true for as 
the extra term only changes the mean. 

Measurements are made on the output beam by split- 
ting the beam into -f 45° and —45° linearly polarized com- 
ponents at a polarizing beam splitter and measuring the 
intensities at the two outputs. We then define the pho- 
tocurrent as 

Ic{t) = lim , (13) 



^5t 



where and are the photon counts at the +45° 
and —45° linearly polarized outputs and v is the photon 
flux Sn/St. We may interpret Idt) as the measured value 
of X divided by VSi, so the photocurrent is given by 

(14) 



/,(^) = a^^(Jf)), + e(^), 



where the subscript c indicates the average for the condi- 
tioned evolution, and ^(t) is a real Gaussian white noise 
term such that (C(i)C(i')> = - 

We take the limit of small time intervals St, though 
for finite photon flux this would mean that the 
SU(2)^IIW(2) contraction is no longer valid, as Sn 
would go to zero. In order for the limit to be rigorously 
correct, we need to take the limit of large ly as well as 
small St. Nevertheless, this limit should be an accurate 
approximation for large photon flux. The conditioned 
master equation is then 

dp,, = TV[4^^^]pcdt + VTdW{t)n[4^^]pc, (15) 

where F — a^v/A, dW{t) = ^{t)dt is an infinitesimal 
Wiener increment, 'D[r]p — rpr'' — {r^rp + prV)/2, and 
'H[r]p = rp + pr'^ — Tr[(r + r'^)p\p. 



III. FEEDBACK 

The quadrature ([l^ ) is measured by detecting the pho- 
tocurrent (p^), thereby yielding TMSS states. The draw- 
back is that this will yield nonzero values of (ji^'') and 
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{Jy~^), whereas for the optimal TMSS states the 
expectation values are zero. In order to obtain states 
closer to the optimal TMSS states, we apply feedback 
to bring these expectation values closer to zero. For ex- 
ample, when cos(f2t) = 1, we may apply a Hamiltonian 



proportional to Jy^'' in order to bring {Jz^') towards 
zero, analogous to the case considered by TMW. When 
sin(rit) = 1, we wish to apply a Hamiltonian propor- 



(+)\ 



tional to ji •* in order to bring {Jy ') towards zero. For 
arbitrary t, we will apply the Hamiltonian 

iffb - F{t)h{t). (16) 

where 

FW^^J^W, (17) 



(-)\ 



and 



) {t) = 4+) cos(m) - sin(m). (18) 



We choose this definition for Jy^^ because [ J- 



— iJx^^ [we will omit the {t) from this point for brevity]. 
As {Jx^'') is close to 2j for the states we are consider- 
ing, this means that rotation about Jy^^ will alter the 
expectation value of ji^^ . This Hamiltonian may be im- 
plemented by a radio frequency magnetic field |Q . This 
generates a Hamiltonian proportional to Jy, which be- 
comes J^^"* due to the Larmor precession of the atomic 
spin. 

Using this feedback, we obtain the unconditioned mas- 
ter equation 

p = -t[^{c^F{t) + F{t)c}, p] + V[c - iF{t)]p, (19) 
with 

c^Vfjf). (20) 

In the single-mode case TMW pointed out that the term 
c' F{t)+F{t)c is proportional to JzJy + JyJz, which is the 
two-axis countertwisting Hamiltonian [ [L6[ . In a similar 
way we can derive a Hamiltonian for the two-mode case 
that produces TMSS states and is analogous to the two- 
axis countertwisting Hamiltonian. It is straightforward 
to show that 

c^Fit) + F{t)c = A(t){[(j(i))2 + (Jf )2 - {4}^f 



, Oy 



sm 



(2m) + (jWjW+jWjW 



+ Ji2)jf)+ jf)jp))cos(2m) 

+ 2(jWjf +jWjf))}. (21) 

In the limit of large the contribution to the evolu- 
tion from the terms containing sm{2Q,t) and cos(2rii) is 
negligible. Omitting these terms yields 

c^F{t) + F{t)c « 2A(t)(jii) -t- JP), (22) 

which indicates that it is possible to produce two-mode 
spin squeezing using the Hamiltonian 

iJcx jWjf +JWJP. (23) 



A. Simple feedback 

Now we will consider the appropriate feedback 
strength, \{t). It is easily shown that the variation in 

{Jz^^^)c due to the conditioned evolution ([l5| ) is 

Tr(jf )dp,) = 2VfdW^[((ji"))2>, - (Jf (24) 
The variation in (ji^^^)c due to the feedback is 



(25) 



If {ji^^)c = 0, then the appropriate feedback to keep this 
equal to zero is 



2r((jfV), 

{'Jx 



(26) 



If the unconditioned state is close to being pure, then 
we may also use this expression with the unconditioned 
averages without significant loss of accuracy: 



m 



2r((jf))2) 



(+)\ 



(27) 



Note that this result is comparable to that found for 
the case of feedback for single-mode spin squeezing by 
TMW: 



2r(j| 



(28) 



This is due to the fact that the commutation relations 
for Ja;, Jy, and Jz are identical to those for Jx'^\ Jy^\ 
and J; 



(O) 



B. Analytic feedback 

Next we will consider the possibilities for using an an- 
alytic expression for the feedback strength A(t). In Ref. 
B the approximate relations 



(j|)«(4R + 2/jr\ 

(J,> «jexp(-ri/2), 
are used to derive the analytic feedback scheme 

rexp(rt/2) 



m 



1 + 2jTt 



(29) 
(30) 



(31) 



In the two-mode case, it is possible to obtain the corre- 
sponding relation for (ji'*'^), but not for ((ji^-*)^). 



4 



In these derivations it is convenient to define the scaled 
variables v = Tt and A(t) = X{t) /T. In terms of these 
variables the master equation (n9|) becomes 



iA{t) 



(32) 



where the dot indicates a derivative with respect to v. 
Using these variables, the evolution is independent of F. 
The value of T does not qualitatively affect the evolution 
and merely provides a scaling to the time and to X(t). 
It is straightforward to show that 

|^(4+))=Tr(j(+)p) 

= -i(4+))+2A((jf))^)-iA2(4+)) 
2((jf))2)2 



(33) 



Initially this derivative is equal to zero. However, the 



value of ((4^"*)^) falls rapidly and that of {Jx'^') only 
falls slowly, so the second term quickly becomes negligi- 
ble. If we neglect the second term, the evolution of (4"*^') 
is 



(+)\ 



(4+))«2jexp(-«/2), 



(34) 



which is equivalent to that in the single-mode case. This 
technique fails to find the corresponding expression for 
((4^'')^), but we have found numerically that a good 
approximation is given by 



((4^))^ 



-v/4 



2v + l/j' 



(35) 



Using these results in the expression for the feedback 
(p7|), we obtain 



rexp(R/4) 
l + 2jTt ' 



(36) 



C. Optimal feedback 

The third alternative that we will consider is optimal 
feedback. In order to derive this, it is convenient to define 
the variables 



c 



((4+')^ + (4''n 

2j 



(J 



(+)\ 



X 



2j 



(37) 



The initial state used is that with individual coherent 
spin states in each arm, so both C and x equal to 
1. As time progresses both x s-^d C decrease. In order 
for the feedback to be optimal, the value of ( should be 
as small as possible for each value of x- Therefore, if x 
decreases by an amount Sx, then the amount by which 



C decreases should be maximized. This means that the 
feedback that is optimum will be that which produces 
the maximum slope. Therefore, in order to determine 
the optimal feedback, we wish to solve 



dAdx 



This can be solved using 



d 

dA 



(+)n2 



(4"')^> 



(+)\ 



= 0. 



(38) 



(39) 



Using the master equation 
to show that 



9|), it is straightforward 



Tr[(4+))p] =-i[l + A2](4+))+2A((4"))2). (40) 

Determining the corresponding expression for the numer- 
ator in Eq. (|3^) is more difficult, and if it is performed 
directly it leads to an extremely complicated expression. 
However, there is a more convenient way of determining 
this expression. Note that 



sm' 



(41) 



-I- 4+^4"') cos(rJi) sin{nt)). 
For large we may average over which gives 

((4"'n-i((4+'f + (4"'n. (42) 

Numerically this approximation was found to be very ac- 
curate. Using this result we find that 

iTr[((4+))^ + (4-))2)p]«Tr[(4^^))^p] 

= A^ji+^f iJi^^yf) - A(44^^)4+'4"' + 4+^ 

(43) 

Thus we find that Eq. (|3^) becomes 
A^d - Ae 



_d_ 
dA 



-/ - AV + Ag 



= 0, 



(44) 



where 



d = ((4+'f-(4^^'f>, 
e = (44")4+)4^)-H4+^ 
/-^(4+4 

5 = 2((4"))2). 
Solving this gives 



A 



-fd±^ifd)^ + ef{fe~dg) 



(45) 



(46) 



fe - dg 

Numerically it is found that the correct solution that 
maximizes the slope is that with the positive sign. 

Note that this derivation of the optimum feedback does 
not rely on any commutation relations other than the 
usual su(2) commutation relations. It therefore is also 
applicable to the case of feedback for single-mode spin 



squeezing, with J; 
and Jj. 



(+) T(n) 



4"^ and 4^"* replaced by J^, J-, 
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FIG. 1: The values of C, (continuous line), x (dash-dotted 
Une), and the purity (dashed line) as a function of time for 
simple feedback (^7|) with a spin of j = 5. The vertical dot- 
ted line indicates the time at which the state is close to the 
optimal state considered for teleportation in Ref. H. 



IV. NUMERICAL RESULTS 

A. Unconditioned master equation 

The unconditioned master equation was solved 
using a simple, finite step method with step sizes of 
5v = 1/1000. It was found that better results were ob- 
tained as Q was increased. In order to estimate the best 
result in the limit of large Q, was assigned the maxi- 
mum value possible with this step size, 'k/2Sv. The initial 
condition used was that where the individual modes were 
in independent coherent spin states oriented along the x 
axis. 

At each time step the quantities x, and the purity 
P = Tr[p2] (47) 

were calculated. The results for the simple feedback 
scheme of Eq. (|2^) and a spin of j = 5 are shown in 
Fig. Initially the system behaves as we would expect, 
with the sum of the variances decreasing, until a time 
of approximately u = 3. Then the variances dramatically 
increase, and the purity drops almost to zero. After this, 
however, the system stabilizes. Note also that the value 
of X drops regularly as time increases until the change at 
w = 3. 

Recall that the system must be entangled if 

((Ji+^)^ + {jf?) < (J^^'^'') §. which is equivalent to 
C < X- Therefore a state with C < X can be described 
as TMSS. Evidently the feedback produces quite strong 
spin squeezing. In order to determine how close to opti- 
mal the states produced by this feedback are, the value 
of C was plotted against x ^-nd compared with the plot 
for optimal TMSS states in Fig. || As time progresses 



FIG. 2: The values of ( plotted against x for spin j = 5. The 
relation for optimal states is shown as the continuous line, the 
results for the simple feedback of Eq. (hi t} as the dash-dotted 
line, the analytic feedback of Eq. (p6p as the dashed line, 
optimal feedback as the dotted line, and the countertwisting 
Hamiltonian ( [23| ) as the pluses. The cross shows the position 
of the optimal state considered for teleportation in Ref. M . 



the state travels from the upper right corner towards the 
lower left corner of the figure. 

The state starts out very close to optimum and does 
not deviate significantly from optimum until < 0.4. 
Eventually the value of C increases dramatically, and the 
system spirals towards its equilibrium state at (0.5,0.5). 
The results for j — 5 are typical, and similar results are 
obtained for other spins j > 1/2. 

For the purposes of teleportation, the important issue 
is how close it is possible to get to the optimal TMSS 
states considered for teleportation in Ref. [|| . This state 
is indicated by the cross in Fig. As can be seen, the 
states obtained by feedback closely approach this optimal 
state. Similar results are obtained for higher spin, but for 
smaller spin the states obtained by feedback are further 
from the TMSS states considered for teleportation. 

The plots of C versus x for analytic feedback and opti- 
mal feedback are also shown in Fig. ^. The results for an- 
alytic feedback are very similar to those for simple feed- 
back, with dramatically increasing at later times. In 
contrast, under optimal feedback the value of C mono- 
tonically decreases to an asymptotic value for x = 0. For 
small times the three feedback schemes produce very sim- 
ilar results, and the states obtained using optimal feed- 
back do not pass significantly closer to the optimal state 
for teleportation. 

The last alternative that we consider in this section is 
the Hamiltonian (^3|) , which is analogous to the two-axis 
countertwisting Hamiltonian in the single-mode case. 
The results for this Hamiltonian are also shown in Fig. ||. 
For this case there is significant spin squeezing, but the 
squeezing is generally poorer than for feedback. Never- 
theless, at later times the value of ( does not rise dramat- 
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FIG. 3: The values of plotted against x for spin j ~ 5. The 
relation for optimal states is shown as the dashed line, the 
results for the feedback of Eq. (^6|) are shown as the continu- 
ous line, and without feedback as the dotted line. All results 
are for the conditioned master equation. The cross shows the 
position of the optimal state considered for teleportation in 
Ref. |]. 



cursions from it. In particular, the value of C does not 
increase dramatically at large times. This indicates that 
improved results may be obtained using feedback based 
on the conditioned state. 

On the other hand, implementing this feedback may 
be far more difficult. In the forms of feedback considered 
for the unconditioned equation, the feedback was propor- 
tional to the measurement record, with a proportionality 
constant A that is a function of time only. This feedback 
may be performed efficiently by calculating X{t) before- 
hand and programming it into a field programmable gate 
array In contrast, the value of A given by Eq. ( |26| ) 
is dependent on the measurement record, and therefore 
would need to be calculated during the measurement, in- 
troducing a significant time delay. 

The results without feedback quickly diverge from the 
line for optimal states, although they do not differ from 
it greatly. The major difference is at later times. While, 
with feedback, the value of C, is reduced very close to 
zero, without feedback the value of ^ cannot be reduced 
as far. This means that, although a significant reduction 
in the variances can be achieved via feedback for strong 
QND measurements (where Ft > 1), feedback does not 
give much improvement for weak measurements. 



ically (until x falls below zero) , as opposed to the results 
for the simple or analytic feedback. 



V. COMPARISON WITH SINGLE-MODE CASE 



B. Conditioned master equation 

Next we consider the results for the conditioned master 
equation. There are two examples where it is necessary to 
consider the conditioned master equation. In the case of 
simple feedback, at the time that increases dramatically 
the purity of the state has greatly decreased. This means 
that the feedback based on the unconditioned averages, 
Eq. (p7|), is a poor approximation of the feedback based 
on conditioned averages, Eq. (|2^), and will not accurately 

keep the average ( ji^'')c equal to zero. For this reason we 
may obtain better results using the conditioned master 
equation and the feedback (^6|). 

The second case that we consider here is that without 
any feedback. In this case there is no reduction in the 
variances for the unconditioned equation. If we consider 
the conditioned equation, however, there is a reduction in 



the variances, though the means {Jz )c and {Jy )c are 
nonzero. For this reason wc will consider the conditioned 
master equation for this case, with ^ defined by 

C- ^- . (48) 

To minimize the random differences between the two 
cases the same random numbers were used for each. The 
results for these two cases are plotted in Fig. |^. 

The results for the feedback (^6|) are extremely close 
to the line for optimal states, with only temporary ex- 



There is clearly a great similarity between the single- 
mode and two-mode cases. The operators ji"'''' , Jy^"^ , and 

ji^^ obey exactly the same commutation relations as Jx , 
Jy, and Jz in the single-mode case. This means that, for 
example, the optimal feedback also applies to the single- 
mode case (as was mentioned in Sec. Ill C| ) . As we may 



replace the sum of the variances {{Jz)'^ + {Jy 'Y) with 
2{{jf^^Y) , it might appear that this case is identical to 
the single-mode case. Nevertheless, there is a subtle dif- 
ference due to the fact that the operator ji^^ is time 
dependent. 

For the case of spin 1/2, there is complete equivalence 
between the single- mode case (for spin j = 1) and the 
two-mode case. We will take the Hilbert space for the 
single-mode case to be that for two spin 1/2 particles, 
so that it is the same as for the two-mode case. The 
optimal TMSS states for j = 1/2 are also optimal single- 
mode squeezed states in terms of the total spin. That is, 
they minimize {{ji^^Y) for a given {J^^). The converse 
is also true: the states optimized for single-mode spin 
squeezing are automatically optimized for two-mode spin 
squeezing. 

In addition, both the single- and two-mode counter- 
twisting Hamiltonians produce optimal spin squeezed 
states. Not only this, but the feedback as given by Eq. 
( p8| ) in the single-mode case, or Eq. ( p7| ) in the two-mode 
case, and the optimal feedback described in Sec. Ill C| , 
give optimal spin squeezed states. The only feedback 
that does not give optimal states is the analytic feedback 
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FIG. 4: The values of ( plotted against x for spin j — 1/2 in 
the two-mode case and j = 1 in the single-mode case. The 
relation for optimal states, single- or two-mode countertwist- 
i ng. the feedbacks of Eq. (ESI) in the single-mode case and Eq. 
(P7|) in the two-mode case, and the optimal feedback of Sec. 
[Ill Q are shown as the continuous line. The results for the 
analytic feedback of Eq. ( p6| ) in the two-mode case, and Eq. 
( ^H ) in the single-mode case are shown as the dashed line and 
dash-dotted line, respectively. 

for the single- or two-mode case. These results are de- 
picted in Fig. IJ. In this figure the variables C and x are 
defined for the single-mode case as 

C = 2(j2)/j- x=(J.>/j. (49) 

To explain these results analytically, first note that, for 
spin 1/2, 

JW 4+) + 4+) Ji+) = 2(j(i) + J^jf^). (50) 

This means that the count ert wist ing Hamiltonians for 
the one- and two-mode cases are identical, and therefore 
produce identical states. 

In addition we find that (ji"*"^)^ — {Jy ^)^ commutes 
with these Hamiltonians. This means that for the states 
produced by countertwisting, the values of {{Jz^^Y) and 
((4"^)2) win be identical. This means that, if the states 
produced are optimal single-mode spin squeezed states 
(i.e., {{Jz^^Y) is minimized), then they must also be op- 
timal TMSS states. 

To show that optimal single-mode spin squeezed states 
are produced by the countertwisting Hamiltonian, note 
that for total spin 1 we have the differential equations 

|-(J.) = 2(j2-4), (51) 

|-(jf) = -(J.), (52) 

^{JD = (J.)- (53) 
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X 

FIG. 5: The values of ( for optimal states plotted against x 
for total spin 2 for the two-mode case (continuous line) and 
the single-mode case (dotted line), and for total spin 10 for 
the two-mode case (dashed line) and the single-mode case 
(dash-dotted line). 

The solution of these equations is 

(J.)2=4((j|)-(jf)2). (54) 

This is the relation for optimal single-mode spin squeezed 
states as given by S0rensen and M0lmer [Q . This shows 
that the states produced by the one- and two- mode coun- 
tertwisting Hamiltonians and the one- and two-mode op- 
timal states are all identical. 

For the case of feedback, it is sufficient to show that the 
feedbacks of Eqs. (^8|) and ( p7| ) give optimal states, as this 
will imply that the optimal feedback gives optimal spin 
squeezed states. The derivation in this case is lengthy, 
and is given in Appendix ^ It is interesting to note that 
in this case the optimal feedback is given by Eq. (|B7|), 
and differs dramatically from the analytic feedback used 
for higher spin. 

The equivalence between the single- and two-mode 
cases does not hold for any total spin higher than 1. 
As shown in Fig. ||, there are significant differences be- 
tween the optimal single- and two-mode spin squeezed 
states for total spin above 1. As can be seen, for a total 
spin of 2 or more the values of ( are significantly less for 
the single-mode case than for the two-mode case. This 
reflects the fact that in the two-mode case we wish to 
minimize ((4"'')^) as weh as ((^i'*'^)^). 

As there are such strong similarities between the single- 
mode and two-mode cases, it is reasonable to apply the 
squeezing parameter as considered in Refs. |l^, |l^ to 
the two-mode case. In the single-mode case the squeezing 
parameter was defined by 

= 2j(J.^)/(J.)^. (55) 

With the definitions of ( and x for the single-mode case 
above, this can be expressed as = C/x^- It would seem 
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FIG. 6: The minimum squeezing parameter ^min multiplied 
j + i as a function of j in the two-mode case (a) and 
the single-mode case (h) for optimal states (continuous line), 
the feedbacks of Eqs. (|2^) and ( psj ) for (a) and (b), respec- 
tively (crosses), analytic feedback (pluses), optimal feedback 
(circles), and countertwisting (dotted line). 



reasonable to use the same definition for the two-mode 
case. 

For the two-mode case we find that the minimum value 
of namely C^jn, varies with spin j as shown in Fig. 
^(a). In this figure the values of have been multi- 
pUed hy j + 1 in order to more easily compare the results 
for different spins. Very similar results are obtained for 
the simple feedback of Eq. (|7|) and the optimal feedback 
of Eq. (|4^). Using the analytic feedback of Eq. (^) gives 
results that are similar to, but slightly above, those for 
the other two feedback schemes. All of these feedback 
schemes give results that are significantly above the re- 
sult for the optimal states. Using the countertwisting 
Hamiltonian (g3|) gives results that are higher than those 
for feedback. 

Similar results are obtained for the single-mode case 
[see Fig. ^(b)]. The results for feedback are notice- 
ably above those for optimal states, and the results for 



two-axis countertwisting are significantly above those for 
feedback. To summarize, the scaling constants for each of 
the cases are as given in Table |. In this table we can see 
the similarities between the single-mode and two-mode 
cases: the scaling constants obtained by feedback are sig- 
nificantly above those for optimal states, and the scaling 
constants for two-axis countertwisting (or its two-mode 
equivalent) are significantly above those for feedback. In 
both cases there are only small differences between the 
results for the different types of feedback. The only qual- 
itative difference between the scaling constants for the 
single- and two-mode cases are that in the single-mode 
case the analytic feedback gives slightly better results 
than the simple feedback of Eq. (p7|). In contrast, in the 
two-mode case the scaling constants are indistinguish- 
able. 



VI. EXPERIMENTAL PROSPECTS 

There are some complications in applying this the- 
ory to the experiment of Ref. In this experiment, 
a « 5 X 10"", and w 2 X lO^s"^ This means that 
r « 1.4 X 10~^s~^, which implies that the time required 
for maximal squeezing is on the order of lO^s (over 20 
yr). Therefore a far more intense beam or a stronger in- 
teraction would be required to obtain the two-mode spin 
squeezing described here. 

It must be emphasized that this long time is required 
for the measurements, not the feedback. Feedback may 
be applied to the states obtained in Ref. |^ with a negligi- 
ble increase in the time required for the experiment. Nev- 
ertheless, for the conditions of this experiment, although 
the feedback will bring the means of Ji+^ and 4+^ to- 
wards zero, it will not significantly reduce the variances. 
It is only when the QND measurement can be performed 
strongly enough (i.e., with a strong enough interaction 
over a long enough time) to obtain spin squeezing close 
to maximum that an improvement in the variance is ob- 
tained by using feedback. 

Another problem is spontaneous losses due to absorp- 
tion of QND probe light. The loss rate due to this is 
iV7(7^n/4A^, where TV = 4j is the total number of atoms 
and g is the one-photon Rabi frequency. As in the case of 
single-mode spin squeezing Q , this loss will be very large 
over the time period l/F for free space. Both this prob- 
lem and the problem of the long interaction time may 
be overcome by performing the experiment in a cavity in 
the strong-coupling regime. 

Two other common experimental problems are ineffi- 
cient detectors and time delays. As for the case of single- 
mode spin squeezing, inefficient detectors do not affect 
the scaling. The system has the potential to be far more 
sensitive to time delays than in the single-mode case, 
however. The problem is that, as the spin component 
that is being measured is rotating rapidly, the feedback 
may be correcting for measurements of a different spin 
component. 
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TABLE I: The scaling constants for the minimum squeezing parameter for both single- and two-mode spin squeezing, 
for optimal states, simple feedback, analytic feedback, optimal feedback, and two-axis countertwisting (or its equivalent in the 
two- mode case). 





Optimal 


Optimal feedback 


Analytic feedback 


Simple feedback 


Countertwisting 


Single mode 


1 


1.6468 


1.6593 


1.6777 


1.9562 


Two mode 


3/4 


1.0584 


1.0692 


1.0692 


1.292 



In order to correct the right spin component, the ro- 
tation of the spin component that is measured should 
have completed an integral number of rotations during 
the time delay. Experimentally, this means that the mag- 
netic field should be adjusted such that the frequency 
is an integral multiple of iir/r, where r is the time de- 
lay. Provided that this is done, time delays should not 
be more of a problem than in the single-mode case. 
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VII. CONCLUSIONS 

Two-mode spin squeezed states are important states 
to produce for quantum teleportation. Here we have 
shown that it is possible to produce states very close to 
the optimal TMSS states by adapting the feedback for 
single-mode spin squeezing considered by TMW. These 
states are not conditioned on the measurement record, 
in contrast to the conditional two-mode spin squeezing 
discussed in Ref. . 

Using the simple feedback scheme ( p7| ) it is possible to 
obtain states that are quite close to optimal TMSS states 
for small times, but strongly diverge from these states at 
later times. In particular, for spins above about 5 it is 
possible to obtain states very close to the TMSS states 
considered for teleportation in Ref. |^ . An analytic feed- 
back scheme (^6|) also gives similar results. This feedback 
scheme is more practical experimentally, as the appropri- 
ate feedback strength to be used is easily calculated. 

We have derived the optimal feedback that produces 
the maximum possible spin squeezing. This feedback is 
also applicable to the case of feedback for single-mode 
spin squeezing. This feedback gives states even closer to 
the optimal TMSS states. 

We have also derived a Hamiltonian for the two-mode 
case that is equivalent to the two-axis countertwisting 
Hamiltonian introduced in Ref. |16 . This Hamiltonian 



produces spin squeezing, but not as strongly as the mea- 
surements with feedback. 

In the case of spin 1/2 both the feedback (except for 
the analytic feedback) and the countertwisting Hamil- 
tonian produce optimal TMSS states. These states are 
equivalent to optimal single-mode spin squeezed states if 
the two modes are considered as a single spin-1 system. 
In the single-mode case optimal spin squeezed states are 
produced both by feedback and by the countertwisting 
Hamiltonian. 



APPENDIX A: RELATION OF VARIABLES TO 
EXPERIMENT 

Here we give further explanation of the variables in- 
troduced in Eq. (|l|) . The coupling constant a is given by 
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A{I + 1/2) A 



(Al) 



where a is the resonant absorption cross section, A is the 
area of the transverse cross section of the light beam, / 
is the nuclear spin, 7 is the spontaneous emission rate 
of the upper atomic level, A is the detuning, and is 
the dynamic vector polarizability. We have omitted the 
bounds on the integral (^ for greater generality, so we 
may apply this expression to multiple samples. Explicit 
bounds are unnecessary, as there is no contribution to 
the integral from regions where there are no atoms. 

The continuous spin operators for the sample are de- 
fined as 



Jfc(z,<) 



lim — 

Sz^o 6z 



EX. 



(A2) 



where k € {x, y, z} and is the Pauli operator for parti- 
cle /i. The sum is over all particles between z and z + Sz 
over the cross section of the sample. The operator for 
the entire sample is obtained by integrating over z. 

The field is described by the continuous-mode annihi- 
lation operators ak{z^t), where k = x and y for the x 
polarized and y polarized modes, respectively. The in- 
stantaneous Stokes parameters are 

Sx{z,t) = i[aj,(z,t)a^(z, i) - al{z,t)ay{z,t)], 
Sy{z,t) = i[aj,(z,t)ay(z,i) + a)y{z,t)aj:{z,t)], 

Sz{z,t) = -^[al{z,t)ay{z,t) - al{z,t)a^{z,t)]. (A3) 

The Stokes vector for the entire pulse at position z is 
given by 8(2) = JS{z,t)dt. 
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APPENDIX B: FEEDBACK FOR TOTAL SPIN 1 

Here we show that the feedback of either Eqs. ( p8| ) or 
( p7| ) gives optimal spin squeezed states for a total spin of 
1. To see this, note first that the first term in Eq. ( p^ ) 
is just the same as that produced by the countertwisting 
Hamiltonian, and so will produce optimal states. In order 
to show that the feedback produces optimal states, it 
remains to be shown that the second term, I?[c — iF]p, 
does not alter the evolution. Specifically, in the single- 
mode case 

Tr {,hV[J, - iAJy]p} = -i(l + A2)(J,> + A(J2 + j^), 

(Bl) 

If the state p is an optimal state, then (J^ + Jy) 
this simplifies to 



1, so 



Tr { J,I?[J, - iAJy]p} = -i(l + A2)( J,) + A. (B2) 
Similarly we can show that 

Tr { Jfl?[j, - tAJy]p} - -|(J,) + A2(j2 - J2)_ (B3) 

For optimal states, (Jj) = 1, so this becomes 



Ti {J^,V[J, - iAJy]p} 



A 



(J.)+A2(l-(jf». (B4) 



It is simple to show that both Eqs. ( ]B2| ) and (B4) will be 
zero if Eq. (54) is satisfied, and the feedback is given by 

Therefore, if the state is already in an optimal state, 
and the feedback as given by Eq. is used, then the 



state will continue to be in an optimal state. On the other 
hand, if some other feedback is used, then optimal states 
will not be obtained. For example, the analytic feedback 
considered by LMW does not give optimal states. In or- 
der to determine analytic feedback that will give optimal 
states, note that the differential equations for (Jx) and 
(J2) are 



dv 
d 

dv 



Jt) = 



- A(2( J2) - 1), 
A 



(B5) 



Solving these using the feedback (M) and Eq. (M) gives 



{Jl) 



-2v 



2^ 



(B6) 



This implies that the value of A should change with time 
as 



Aiv) 



1 



(B7) 



This is quite different from the analytic expression ap- 
plied for larger spins. 

The case for feedback for two-mode spin squeezing is 
analogous to this, except that Jx, Jy, and Jz are replaced 

with ji^\ Jy^\ and ji^^ This therefore shows that 
optimal states are obtained in the one- and two-mode 
cases using the feedback given by Eqs. (^) and (|27" 
respectively. 
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